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ABSTRACT 

We present results from observations of the most famous starburst galaxy M82 
with the High-Resolution Camera onboard the Chandra X-Ray Observatory. We 
found nine sources in the central 1' x 1' region, but no source was detected at 
the galactic center. Comparing the observations on 1999 October 28 and on 2000 
January 20, we found four of the nine sources showed significant time variability. 
In particular, CXOM82 J095550. 2+694047, which is 9" away from the galactic 
center, showed extremely large time variability. We conclude that this source 
is the origin of the hard X-ray time variability of M82 detected with ASCA. 
Assuming a spectral shape obtained by the ASCA observation, its luminosity in 
the 0.5 - 10 keV band changed from 1.2x10*^*^ erg s^^ on 1999 October 28 to 
8.7x10^° erg s~^ on 2000 January 20. 

Subject headings: galaxies: starburst — galaxies: active — galaxies: individual 
(M82) — X-rays: galaxies 

1. Introduction 

The X-ray spectrum of the starburst galaxy M82 measured with ASCA (Tanaka, Inoue, 
& Holt 1994) consists of soft, medium, and hard components (Tsuru et al. 1997). The hard 
component, which dominates the X-ray spectrum above 2 keV, was found to show time 
variability by monitoring M82 with ASCA in 1996 (Matsumoto & Tsuru 1999; Ptak & 
Griffiths 1999). Assuming a distance of 3.9 Mpc (i.e. 1" - 19 pc) (Sakai & Madore 1999), 
its luminosity in the 0.5 - 10 keV band changed between 4.5x10^° erg s^^ and 1.6x10^^ 
erg s~^ at various time scales from 1x10"^ s to a month. The most plausible explanation for 
the variability is that a low-luminosity AGN exists in M82. If this explanation is correct. 
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M82 would be one of the most important objects for the study of a relation between an AGN 
and starburst activity (e.g. Umemura, Fukue, & Mineshige 1997) because of its proximity. 
However, hard X-ray spectra of M82 obtained with Ginga and BeppoSAX, which can detect 
the X-ray photons of higher energy than ASCA, differ from typical AGN spectra, since they 
cannot be fitted by a power-law model (Tsuru 1992; Cappi et al. 1999). Therefore, it is also 
possible that the origin of the hard component is a new type of compact X-ray source. 

Furthermore, Matsumoto & Tsuru (1999) compared the ASCA image with the ROSAT 
HRI image (Strickland, Ponman, & Stevens 1997), and showed that the hard component 
probably comes from the X-ray brightest source detected with the ROSAT HRI. The HRI 
source was found to show time variability (Collula et al. 1994). Since Stevens, Strickland, & 
Wills (1999) showed that the source is away from the dynamical center of M82, this galaxy 
may harbor an "off-center" AGN. 

Thus, the source of the hard component of M82 would be a quite interesting and 
important object. For further investigation, it is necessary to determine the position of 
the hard component precisely and to find a counterpart in other wavelengths if possible. 
Therefore, we analyzed the data of M82 obtained with the High-Resolution Camera (HRC) 
(e.g. Murray et al. 1987, 1997) onboard the Chandra X-Ray Observatory ('CX(9j(Weisskopf 
et al. 1995). The HRC has sensitivity in the 0.1 - 10 keV band with a peak at about 1 keV, 
and it has the highest angular resolution (FWHM ~ O'.'S) of all instruments onboard X-ray 
observatories to date. Uncertainties in this paper refer to la confidence limits. 

2. Data Analysis and Results 

M82 was observed twice with CXO as a calibration target using the HRC-1 on 1999 
October 28 and 2000 January 20. The HRC data were processed with the standard 
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procedure by the Chandra X-ray Center (CXC) ^°, and we used only the fully processed 
science products (level = 2) as event files. We also applied an event screening procedure 
developed by the HRC team to eliminate "ghost" events. The exposure times are 2788 s for 
the first observation and 17684 s for the second. 

We found that all X-ray bright sources which can be the origin of the ASCA hard 
component exist within the central 1' x V (= 1.1 kpc x 1.1 kpc) region, which is shown in 
Figure 1. The position of the bright off-center source found with the ROSAT HRI (Stevens, 
Strickland, & Wills 1999) is also within the field. Therefore, we concentrate on the central 
V X 1' region. In Figure 1, we see nine prominent sources, which are designated with 
circles and numbers. More detailed analysis including diffuse emission and other fainter 
sources will be presented in a forthcoming paper (M. Ward et al., in preparation). We 
determined the positions of the peaks of the nine HRC sources with the wavelet algorithm 
using Chandra Interactive Analysis of Observations (CIAO), and named the HRC sources 
by using these positions (Table 1). The position uncertainty is less than ~ O'.'l. 

According to "the CXC memo on astrometry problems" ^^, the HRC event files we 
used could have offsets up to 10" in the celestial coordinates. To check the reliability of the 
coordinates in the event files, we compared the coordinates of the HRC sources with the 
Two Micron All Sky Survey (2MASS) Point Source Catalog. There are 17 2MASS sources 
in the field of Figure 1. According to the CXC, the offset of the celestial coordinate in 
ordinary CXO event files is 1" (RMS). The position uncertainty in the 2MASS Catalog is 
~ O'.'l. Considering these uncertainties, we found that three 2MASS sources out of 17 agree 
with the HRC sources (Nos. 1, 6, and 8). Since the probability of a chance coincidence is 
less than 0.1 %, we can assume that the HRC coordinates are reliable. We then compared 



^°http://asc. harvard.edu/ 

^^http://asc. harvard.edu/ciao/caveats/aspect4. html 
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the coordinates of the HRC sources with the 5 GHz sources (Muxlow et al. 1994) whose 
position accuracy is ~ C'S. The fact that four HRC sources (Nos. 4, 5, 6, and 7) have 
counterparts on the 5 GHz map also supports the rehabihty of the HRC coordinates. The 
identifications of the HRC sources are shown in Table 1. 

We estimated the counting rates of the HRC sources using the X-ray events within 
the circular regions shown in Figure 1. The backgrounds including the diffuse emission 
were estimated using source-free regions around the sources. The counting rates are listed 
in Table 1. Four HRC sources (Nos. 5, 7, 8, and 9) show significant time variability. In 
particular, No. 5 disappeared in the second observation, and No. 7 became brighter by a 
factor of 7. The counting rates of No. 1 and No. 4 did not change significantly between the 
two observations. 

The dynamical center at {a,S)j2ooo = (9'^55™51.9®, 69°40'47.1"), which is determined 
by the radio observation of the motion of H I gas, is shown as the green cross in Figure 1 
(Weliachew, Fomalont, & Greisen 1984). The position error circle with a radius of 2" 
is also shown as the green circle. We should note that all nine HRC sources are clearly 
away from the dynamical center. The counting rate of the dynamical center estimated 
with the extraction radius of 2" is (1.79 ± 1.48) xlO"''^ c s"^ for 1999 October 28, and 
(1.12 ± 6.84) xlO^^ c s^^ for 2000 January 20. Thus, we found no significant source at the 
dynamical center. 

We estimated the HRC counting rate of the ASCA hard component (Matsumoto & 
Tsuru 1999) using W3PIMMS v3.0 ^^: the expected counting rate is 0.72 c s~^ in the 
highest state and 9.7x10^^ c s^^ in the lowest. Only the counting rate of No. 7 is consistent 
with the expected counting rate. Furthermore, the position of the ASCA hard component is 
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consistent with that of No. 7 (Matsumoto & Tsuru 1999). Therefore, we can conclude that 
the variabihty of the ASCA hard component is due to No. 7 (CXOM82 J095551.4+694044). 
The separation between No. 7 and the dynamical center is 9" (~ 170 pc). 

3. Discussion 

We found that 41.5+59.7 (Kronberg, Biermann, & Schwab 1985) is a candidate of the 
radio counterpart of No. 7 as well as 41.30+59.6 (Muxlow et al. 1994). 41.5+59.7 is 0'.'76 
away from No. 7, while the separation between 41.30+59.6 and No. 7 is 0'.'96. According to 
the morphology and spectral shape in the 5 GHz radio band, Muxlow et al. (1994) suggested 
that 41.30+59.6 is a young supernova remnant (SNR). 41.5+59.7 show a 100 % drop in the 
radio flux within a year, and its radio decay time scale and spectrum are very similar to 
SN 1983n (Kronberg & Sramek 1985; Kronberg et al. 2000). Therefore, this source may also 
be a SNR. The hard X-ray emission of M82 was found to show the short-term variability 
(Matsumoto & Tsuru 1999; Ptak & Griffiths 1999). Since it is rather difficult to explain 
the short-term variability in terms of a SNR origin, both radio sources may not be a real 
counterpart of No. 7. 

If we assume the spectral shape of No. 7 is an absorbed thermal bremsstrahlung model 
with a temperature of 10 keV and a column density of 10^^ cm~^, which is the typical 
spectral shape of the ASCA hard component (Matsumoto & Tsuru 1999), the peak X-ray 
flux of No. 7 is 6.6x10"^^ erg s~^ cm~^ in the 2-10 keV band. According to the the 
logA^ — logs' relation (Ueda et al. 1999), the probability that a source as bright as No. 7 
exists in the 1' x V field is ~ 0.3 %. Therefore, No. 7 is probably not a background 
AGN. The ASCA spectrum of the variable source obtained by subtracting the spectrum 
of the lowest state from the highest state shows heavy absorption (the column density is 
~ 10^^ cm~^) (Matsumoto & Tsuru 1999). Since the Galactic absorption toward M82 is 
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4x10^° cm~^ (Dickey & Lockman 1990), the variable source is embedded deeply in M82, 
and hence No. 7 is probably not a foreground source, unless the source has extremely large 
intrinsic absorption. 

Assuming the absorbed bremsstrahlung model of 10 keV and 10^^ cm^^, unabsorbed 
X-ray luminosity in the 0.5 - 10 keV band (L^^^^^^^^) is expressed using an HRC counting 
rate (Crrc) as 



^a5-10kcV^^g7^^Q38grgs-l X 



Q 



HRC 



10-3 C S" 



The L^5-iokcV Qf ^Q^ 7 ^g^g estimated to be 1.2x10^° erg s'^ on 1999 October 28 and 
8.7x10^° erg s^^ on 2000 January 20. 

If we assume that No. 7 is a black hole (BH) and that the maximum luminosity does 
not exceed the Eddington luminosity, the mass of the BH must be larger than 700 Mq, 
and No. 7 is not a stellar-mass BH (~ 10 Mq). Since No. 7 is 9" away from the dynamical 
center, the mass of No. 7 must be much smaller than the gravitational mass within 9" 
from the center which is 4x10^ Mq (McLeod et al. 1993), otherwise the dynamical center 
would be shifted from the current position. Therefore, No. 7 is at the low end of the mass 
distribution of super-massive BHs (10^ - 10^ Mq) or a medium-massive BH (10^ - 10^ 
Mq). The possibility of the medium-massive BH is discussed in detail in Matsushita et al. 
(2000) along with the discovery of an expanding molecular superbubble surrounding No. 7 
(Matsushita 2000). Other possibilities such as an X-ray binary source whose jet is strongly 
beamed at us cannot be excluded. Further investigation including other wavelengths is 
strongly encouraged to reveal the true character of No. 7 (CXOM82 J095550. 2-1-694047). 

If we use equation (1) to the sources other than No. 7, their L^^^^^ ^ s are much 
greater than the Eddington luminosity for a 1.4 Mq object. Therefore, these sources may 
also be BHs. 

Though we found no significant source at the dynamical center, it is still possible that 



a faint X-ray source such as Sgr A* exists (Koyama et al. 1996). Assuming an absorbed 
power-low model with a photon index of 1.7 and a column density of 10^^ cm^^, the upper 
hmit of i^^5-iokcV ^£ ^YiQ dynamical center is 5.4x10^^ erg s~^ for 1999 October 28 and 
1.3 xlO^^ erg s^^ for 2000 January 20. 

This paper makes use of data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/ California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. The authors are grateful to Miss 
Deborah Gage for careful review of the manuscript. We also thank Dr. P. Kaaret for 
valuable comments. H. M. and S. M. are supported by the JSPS postdoctoral Fellowships 
for Research Abroad. 
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Table 1. X-ray sources at the central region of M82 detected with the HRC. 



Number 



Source Name'' 



Identification 



Counting Rate (10 '^ c s ^) 



Infrared'^ 



Radio'^ 



R'' 



1999 Oct 28 2000 Jan 20 



1 CXOM82 J095554.7+694101 0955547+694100 

2 CXOM82 J095552. 3+694054 

3 CXOM82 J095551. 6+694036 

4 CXOM82 J095551. 4+694044 

5 CXOM82 J095551. 1+694045 

6 CXOM82 J095550. 7+694044 0955507+694043 

7 CXOM82 J095550.2+694047 

8 CXOM82 J095547.5+694036 0955475+694036 

9 CXOM82 J095546.6+694041 

dynamical center 





2{'0 


2.51+1.08 


2.42+0.40 




2{'0 


6.10+2.00 


3.15+0.67 




2{'0 


4.30+2.43 


8.43+0.95 


42.65+57.8 


1'.'3 


8.04+2.14 


6.27+0.72 


42.21+59.0 


1"4 


20.9+3.1 


0.292+0.505 


41.95+57.5 


1"1 


2.99+1.54 


5.20+0.65 


41.30+59.6 or 41.5+59.7<= 


2"0 


71.1+5.9 


520+5 




2"0 


24.8+3.1 


4.89+0.66 




2"0 


4.63+1.48 


22.8+1.2 




2"0 


1.79+1.48 


0.112+0.684 



^Sources are named according to the Chandra Source Naming Convention (http://asc.harvard.edu/udocs/naming.html). 
For example, CXOM82 J095554.7+694101 is a source at {a, 5)j2ooo=(9''55™54.7=,69°41'01"). 

''The Two Micron All Sky Survey (2MASS) Point Source Catalog (http://www.ipac.caltech.edu/2mass/) 

'^Muxlow ct al. (1994) 

■^Extraction radius used to estimate the counting rate. The radii for Nos. 4, 5, and 6 were limited to below 2", because they 

are close to each other. 

°Kronberg, Biermann, & Schwab (1985) 

''(",<5)j2000 = (9''55™51.9",69°40'47.1") (Wcliachew, Fomalont, & Greisen 1984) 



Note. — Errors for the counting rates are given for Icr confidence. 
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Fig. 1.— Central 1' x 1' (= 1.1 kpc x 1.1 kpc) region of M82: (a) 1999 October 28, (b) 
2000 January 20. The pixel size is 0'.'26 and the image has been smoothed with a Gaussian 
filter with a a of 1 pixel. The color scale shows the total count per pixel. The green cross 
shows the radio kinematic center, and the green circle is its position error circle with a radius 
of 2" (Weliachew, Fomalont, & Greisen 1984). The sources we analyzed in this paper are 
designated with white circles and numbers. The circular regions were used to extract X-ray 
events for the sources. 



